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ABSTRACT 

The radius of neutron stars can in principle be measured via the normalisation of a blackbody fitted 
to the X-ray spectrum during thermonuclear (type-I) X-ray bursts, although few previous studies have 
addressed the reliability of such measurements. Here we examine the apparent radius in a homogeneous 
sample of long, mixed H/He bursts from the low-mass X-ray binaries GS 1826—24 and KS 1731—26. 
The measured blackbody normalisation (proportional to the emitting area) in these bursts is constant 
over a period of up to 60 s in the burst tail, even though the flux (blackbody temperature) decreased 
by a factor of 60-75% (30-40%). The typical rms variation in the mean normalisation from burst to 
burst was 3-5%, although a variation of 17% was found between bursts observed from GS 1826—24 
in two epochs. A comparison of the time-resolved spectroscopic measurements during bursts from 
the two epochs shows that the normalisation evolves consistently through the burst rise and peak, 
but subsequently increases further in the earlier epoch bursts. The elevated normalisation values may 
arise from a change in the anisotropy of the burst emission, or alternatively variations in the spectral 
correction factor, f c , of order 10%. Since burst samples observed from systems other than GS 1826—24 
are more heterogeneous, we expect that systematic uncertainties of at least 10% are likely to apply 
generally to measurements of neutron-star radii, unless the effects described here can be corrected for. 
Subject headings: stars: neutron — X-rays: bursts — X-rays: individual(GS 1826—24) — X-rays: 
individual(KS 1731—26) — techniques: spectroscopic 



1. INTRODUCTION 

Interest has been raised in recent years in the prospects 
of inferring the neutron-star mass and radius from ther- 
monuclear bursts. Such a possibility can provide strin- 
gent constraints on the neu tron-star equation of state , 
which remains uncertain (e.g. lLattimer &: Prakash 2007). 
From combining measurements of the blackbody nor- 
malisation (from time-resolved X-ray spectral fits to 
the burst spectra), Eddington luminosity, and the dis- 
tance, confidence limits on the mass and radius of several 
neutron stars in low-mass X-ray binaries (LMXBs) has 
been estimated (|Ozel et all l2009t fGuver et all l2010al lbl: 
lOzel et al.1 120111 : see also iSteiner et afl I2010I ). Each 
of these quantities is challenging to measure alone, 
but particularly substantial systematic errors are known 
to affect the blackbody normalisation. The normal- 
isation can rise or fall steadily throughout the burst 
tail, depending roughly on the duration of the burst 
(jBhattacharyya et al.l [2010). Furthermore, in some 
sources different bursts provide sig n ificant ly different 
normalisation values. iGiiver et al.l (|2011| ) quantified 
some of these effects at the 3-7% level, although in an 
earlier analysis of the burst source EXO 1745—248 found 
apparent radii measured from two bursts varying by 11% 
(|Ozel et al.ll2009f) . Additionally, analysis of bursts from 
another globular cluster source, 4U 1724—307, found 
blackbody normalisations which varied by a factor of 
w 2 (i.e. a 40% variation in the in ferred radius) be- 
tween short- and long-duration bursts (jSuleimanov et al.1 
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I2011al) . 

Independent of the issues for measurement, are un- 
certainties about precisely how the neutron star's at- 
mosphere affects the emerging radiation. Although 
the burst spectra are typically found ob servationally 
to be consistent with a blackbody (e.g. iSwank et al.l 
Il977t iKuulkers et al.ll2002j ). scattering effects have long 
been understood to distort the spe ctrum sufficiently to 
bias t he measured temperature (e.g. lLondon et al.lll98"H 
1986). This distortion is usually parameterised via a 
spectral distortion factor f c — T\^/T e s where Tbb is the 
measured blackbody (or colour) temperature, and T e s is 
the effective temperature of the atmosphere. Most recent 
work adopt a narrow range of f c = 1.3-1.4 a t burst lumi- 
nosit ies well below the Eddington limit (e.g. Mad ei et al.l 
|2004| ). and often neglect t he variations in f„ that may 
arise during the burst (e.g. iSuleimanov et al J 12011 bp ). A 
more rigorous approach involves fitting the observed vari- 
ation in the blackbody normalisation ( in response to the 
changi ng f c ) as a function of flux (e.g. ISuleimanov et al.l 
l2011a| ). although the predicted model curves cannot yet 
reproduce the range of observed behaviour. 

Samples of bursts accumulated from individual sources 
can be extremely heterogeneous in their properties. From 
low-duty cycle observations featuring gaps due to Earth 
occultations and other conditions, it is usually impos- 
sible to be confident about the burst recurrence time, 
in which case the detailed ignition conditions, fuel com- 
position, and even the completeness of thermonuclear 
burning are also uncertain. Under such conditions, it is 
difficult to disentangle the various systematic influences 
whic h might influence the normalisation measurements 
(e.g. IGiiver et al.ll20lTl ). Here we investigate the intrin- 
sic reproducability of burst normalisation measurements 
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using a uniform, homogeneous sample of bursts, from the 
low-mass X-ray binaries GS 1826-24 and KS 1731-26. 
We use Rossi X-ray Timing Explorer (RXTE) Propor- 
tional Counter Array (PCA) data to test for intrinsic 
systematic effects which might influence burst normali- 
sation measurements beyond any additional effects which 
might arise from variations in the ignit ion conditions and 
fuel composition. In a related paper, iZamfir et al.ll20ill 
(submitted to ApJ; hereafter Zll) used the same data to 
infer the mass and radius of GS 1826—24. 

2. OBSERVATIONS AND ANALYSIS 

Few bursting sources exhibit trains of bursts with 
consistent lightcurves or recurrence times. The best 
known example is GS 1826—24, so far unique for its 
consistently regular burst behaviour, and high degree 
of uniformity between successive burst lightcurves (e.g. 
iGallowav et al.l 120041 ). Comparison of the burst be- 
haviour and lightcurves suggest that the system ac- 
cretes mi xed helium and hy drogen at roughly solar mass 
fraction ()Heger et al.l [20071. We used observations of 
GS 1826—2 4 taken with the P roportional Counter Ar- 
ray (PCA; Uahoda et al.l 119961 ) onboard the Rossi X- 
ray Timing Explorer (RXTE), from t he catalogue of 
burst s detected over the mission lifetime ()Gallowav et al.l 
l2008t hereafter G08). The flux- recurrence time re- 
lations hip for this sample ha s been extensively stud- 
ied bv I Thompson et all (|2008l ). Optical photometry of 
the mass donor suggests an orbital period of 2.25 hr 
(Meshchcryakov et al. 2010). However, several alias 
peaks are present in the periodogram, and it is possi- 
ble one of these (particularly at 2.05 hr) represents the 
true orbital period. 

We performed a search of the G08 sample for addi- 
tional examples of regular, consistent bursts. Recurrence 
time provides the most obvious way to detect regular 
bursts, although for instruments in low-Earth orbit like 
RXTE, recurrence time measurement is confounded by 
regular interruptions due to occultations of the star by 
the Earth. An alternative approach is to test for con- 
sistency of the burst light curve, via commonly-used pa- 
rameters measuring the duration. The ratio r of the 
fluence Eb to the peak flux, as used by G08, provides 
a simple way of comparing light curves. A scatter plot 
of r against Eb for a source with consistent bursts, will 
show strong clustering, and KS 1731—26 provides the 
next best example after GS 1826—24. Extensive ob- 
servations of KS 1731-26 by the BeppoSAX/WFC sug- 
gest that, at times, this system exhibits regular bursts 
mor e frequently than GS 1826—24 (see e.g. Fig. 3 
from lCornelisse et al.ll2003l ). Despite this, little attention 
has been directed at the burst behaviour of this system. 
One reason is that, unlike GS 1826-24, KS 1731-26 
exhibi ts both radius-exp ansion bursts, and burst oscilla- 
tions (jMuno et al.ll2000f ); previous analyses have focussed 
largely on these phenomena. A series of RXTE obser- 
vations of KS 1731—26 were made in 2000 August and 
September, detecting a total of 14 bursts with highly con- 
sistent lightcurves. WFC observations also made during 
this time show that the burst behaviour was quite reg- 
ular, with recurrence times of between 2-3 hr. An ad- 
ditional example, with a lightcurve very similar to those 
observed in 2000, was detected on 1999 Aug 26. Soon af- 
ter the 2000 observations, the source faded to quiescence, 



and has not been active since (jWiinands et al.ll200l . 

No other sources exhibit clustering as tight in the 
rest of the G08 sample; the next best examples are 
4U 1323-62 and 4U 0836-429. 4U 1323-62 is quite 
well-known for it's regular burst behaviour (e.g. G08), 
and in the RXTE sample, many of the bursts were 
long with t = 28 ± 2 s and similar fluences, of Eb = 
(100±10) x 10 -9 erg cm -2 . However, several of the bursts 
were closely followed by weaker events only a few min- 
utes after the first. This behaviour has be en observed in a 
number of systems (e.g. lKeek et a l. 2010), and is thought 
to arise from delayed ignition of material left over from 
the initial burst. Such events are likely not suitable for 
a detailed test of the consistency of blackbody normal- 
isations between bursts. During a period of activity in 
2003-4, 4U 0836-429 exhibited long (t = 22±4 s) bursts 
at roughly similar fluences, distributed approximately as 
a Gaussian with E b = (270 ± 70) x 10" 9 erg cm" 2 . De- 
spite good coverage of the outburst, the bursts were ap- 
parently not strictly periodic in nature, and there may 
also be contributions to the source flux from a nearby 
HMXB pulsar which is within the PCA field of view. 
Thus, we also exclude this source from our study. 

Where not explicitly stated, the data analysis proce- 
dures are as in G08. Time-resolved spectra in the range 
2-60 keV covering the burst duration were extracted on 
intervals as short as 0.25 s during the burst rise and peak, 
with the bin size increasing gradually into the burst tail 
to maintain roughly the same signal-to-noise level. A 
spectrum taken from a 16-s interval prior to the burst 
was adopted as the background. We re-fit the spectra 
over the energy range 2.5-20 keV using the revised PCA 
response matrices, vll.TQand adopted the recommended 
systematic error of 0.5%. The fitting was undertaken us- 
ing XSpec version 12. In order to accommodate spectral 
bins with low count rates, we adopted Churazov weight- 
ing. No correction for instrumental deadtime was applied 
to the spectra. 

We modelled the effects of interstellar absorption, us- 
ing a multiplicative model component (wabs in XSpec), 
with the column dens ity nn frozen at 4 x 10 21 cm~ 2 
(for GS 1826-24, e.g. lin 't Zand et"aHU999l) and 1 .3 x 
10 22 cm" 2 (for KS 1731-26, e.g. ICackett et all200ll . In 
the original analysis carried out by G08, the neutral ab- 
sorption was determined separately for each burst, from 
the mean value obtained for spectral fits carried out with 
the nn value free to vary. This has a negligible effect 
on the fluxes, but can introduce spurious burst-to-burst 
variations in the blackbody normalisation. We also com- 
puted averaged lightcurves of blackbody spectral param- 
eters for subsets of burst s from GS 1826—24, f ollowing 
the procedure adopted by IGallowav et all (|2004f ). 

3. RESULTS 

We first explored a variety of different approaches to 
measuring the blackbody normalisation in the tail of the 
bursts. We used an event from GS 1826-24 on 2000 
July 1 17:16:37 (#12 in G08). Four of five PCUs (0-3) 
were functioning for this burst. We took the approach 
of iteratively determining the maximum extent in the 
tail (beginning from a few seconds after the burst start) 

4 see \pro tectj|ittp: / /www. universe. nasa. gov/xrays/programs /rxte/pca/ 1 
doc/ rmf /pcarmt- 11.7 
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Table 1 

Different approaches for determining the mean blackbody normalisation 



Treatment 


Mean A"bb 
[(km/^iOkpc) 2 ] 


Time range 

w 


kT range 
(keV) 


Flux range 
[10 — 9 erg cm -2 s _1 ] 


xl (DOF) 


averaged, variable bins 
joint fit, variable binsize 
averaged, 0.25-s bins 
joint fit, 0.25-s bins 


103.3 ±0.7 
107.0 ± 1.1 
103.5 ±0.7 
105.5 ±0.9 


3.0-65.25 
11.50-32.25 
3.25-50.75 
4.75-32.00 


2.32-1.55 
2.16-1.79 
2.32-1.46 
2.32-1.72 


28.69-6.09 
24.26-11.80 
28.69-7.05 
28.69-11.18 


1.45 (95) 
1.13 (942) 
1.28 (186) 
1.08 (2460) 



over which a constant fit to the blackbody normalisations 
was an acceptable fit to below 3<r confidence. We refer 
throughout to the best-fit constant value as (Kbb), to 
distinguish from the individual Kbb values for each time- 
resolved spectrum. 

We found a best-fit value of (Kbb) — 103.3 ± 
0.7 (km/diokpc) 2 over the interval 3.0-65.25 s (relative 
to the start over the burst; Fig. [1]). The reduced 
X 2 / n DOF = xt was 1-45, for 95 degrees of freedom. Re- 
markably, this consistency was maintained despite the 
fitted blackbody temperature kT^b falling from 2.32 to 
1.55 keV. 

In order to double-check our result, we carried out a 
similar procedure to determine the maximum extent over 
which the blackbody normalisation was consistent with a 
constant value, but instead of simply averaging the fitted 
normalisations for each time bin, we performed a joint 
fit in XSpec of the time-resolved spectra to an absorbed 
blackbody model, as used for the individual spectra. We 
froze the neutral column density at the same value as 
used for the individual spectral fits, 4 x 10 21 cm -2 . We 
linked the blackbody normalisation for each of the indi- 
vidual spectra but allowed the blackbody temperatures 
to vary. With this treatment, the fitted -Kbb.joint was 
found over a shorter time interval, between 11.5-32.25 s 
following the burst start, and at a 3% higher value of 
107.0 ± 1.1 (Table [1]). Likely, the higher value is attained 
by omitting some of the spectra in the range 40-60 s, 
when there is a noticeable downward trend in the nor- 
malisation (Fig. [1]). We note that the mean of the nor- 
malisation values over the reduced extent of the burst 
permitted by the joint fits, returns a consistent best-fit 
value, of 106.6 ± 1.1. 

The fact that the two treatments returned a different 
extent of the burst tail over which the normalisation was 
constant may be attributed to the degree of goodness of 
fit of the individual spectra. The distribution of \ 2 for 
the individual fits in the time range 3.0-62.5 s follow- 
ing the burst start is shown in Fig. [2] Compared to 
the expected distribution given the number of degrees of 
freedom (22), the observed distribution is skewed signif- 
icantly to higher values of \ 2 ■ This result suggests that 
residual systematic errors (either instrumental, or possi- 
bly from deviations from a blackbody) prevented a joint 
fit to the same confidence level extending over the same 
extent of the burst tail, compared to the simple average. 

We also explored the effect of the time binning ap- 
proach on the results, by repeating the two analyses 
on data which used uniform 0.25 s bins throughout the 
burst. The bin sizes for the G08 data were 0.25-s from 
the start of the burst through to 14.5 s after the start, 
0.5 s to 24.5 s, 1 s to 43.25 s, and 2 s through to 65.25 s. 
We found that the (Kbb) value obtained with uniform 
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Figure 1. Example burst observed by RXTE from GS 1826-24 
on 2000 July 1 17:16:37 UT (#12 in G08), illustrating the con- 
stancy of the blackbody normalisation over a significant extent of 
the burst tail. Shown is the burst flux (top panel), blackbody 
(colour) temperature fcT^b (middle panel), and blackbody normal- 
isation (in units of (km/dxQj.p;.) 2 ; bottom panel). The extent over 
which the constant-normalisation fit is calculated is illustrated by 
the vertical dotted lines (top and middle panels), and the horizon- 
tal lines (lower panel), with the la error range indicated (dashed 
lines). The reduced-x 2 = xi f° r the constant fit is 1.45, for v = 95 
degrees of freedom. 

0.25-s bins was identical to that with bins of variable 
size, although the extent of the constant fit was slightly 
less (extending to 50.75 s instead of 65.25 s; Table [TJ. 
Similarly, the joint fit to the uniform 0.25-s binned spec- 
tra gave a comparable result to the joint fit for variable 
bins, and over roughly the same extent of the burst. The 
X 2 distribution for the 0.25 s-binned data over the range 
3.25-50.75 s after the burst start was similarly discrepant 
from the expected distribution assuming a statistically 
good fit, but at a lower confidence level (K-S statistic of 
0.12, equivalent to 2.5<j). 

The RXTE PCUs are subject to a short (» 10 /is) 
interval of inactivity following the detection of each 
X-ray photon. This "deadtime" reduces the de- 
tected count rate below what is incident on the de- 
tector (by approximately 3% for an incident rate of 



400 count 



pcir 



For GS 1826-24, the bursts 



peak at w 2500 count s 1 PCU , giving a peak dead- 
time rate of ~ 6% (reducing in the burst tail to w 4% 
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Figure 2. Distribution of fit statistic \ 2 for blackbody fits to time- 
resolved spectra from burst #12 from GS 1826-24, on 2000 July 
1 17:16:37, over the time for which the blackbody normalisation 
was determined to be constant (3.0—65.25 s relative to the burst 
start). The expected distribution for statistically acceptable fits 
with 22 degrees of freedom is overplotted {dotted line). The fit \ 2 
show an excess of high values, as also indicated by the cumulative 
probability distributions (inset); a K-S test returns a value of 0.21, 
indicating disagreement at the 3.4<r level. 

after « 60 s). Since the bursts reach approximately the 
same peak flux, the deadtime correction is also approx- 
imately the same, and can be neglected for comparison 
purposes (we similarly neglect any issues related to the 
absolute PCA calibration). However, since the deadtime 
correction varies with time, it might also be expected to 
result in a slightly different evolution of the normalisation 
throughout the burst. To test this possibility, we exam- 
ined the duration of the consant interval for the burst 
analysed in this section (#12 from G08), with and with- 
out the deatime correction. Although we measured an 
increase in the average normalisation with the deadtime 
correction, as expected, we found no change in the extent 
over which the normalisation was found to be constant. 
Thus, we neglect deadtime corrections for the remainder 
of the analysis in this paper. 

3.1. The full burst sample from GS 1826—24 

We analysed 67 bursts observed from GS 1826—24 by 
RXTE between 1997 and 2007. This sample includes 
the set of 58 complete bursts in G08, excluding seven 
events for which the full lightcurve was not observed, or 
for which no appropriate data modes were available for 
the analysis, and one event which occurred during a slew 
and for which the burst flux could not be determined 
correctlj0. We also analysed an additional 9 bursts from 
observations in 2006 August, which were not part of the 
G08 sample. 

We compared the results of the joint fits Abb joint with 



the average of the individual fits (Abb) for the bursts 
from this sample to determine how the results from our 
sensitivity tests translated to a larger sample. Although 
in some cases the joint fit durations were shorter than 
the constant fit duration for the individual normalisa- 
tions the median duration was 94% of the average in- 
terval. Furthermore, we found no systematic difference 
between the normalisations determined by the two meth- 
ods, and the agreement was at better than 2% for 76% of 
the measurements (maximum deviation was 8.5%). This 
compares favourably to the typical la statistical uncer- 
tainty on the measurements, of 1%. Thus, we conclude 
that the parameter averages (kbb) provide an unbiased 
measure of the blackbody normalisation in the tail of 
these bursts, and adopt that method for further analy- 
sis. 

The (A'bb) for the bursts from GS 1826-24 determined 
from the individual spectral fit parameters (using the 
spectra with variable binsizes; see Q is shown in Fig. [3J 
The typical span of the constant fit was from < 3 s after 
the burst start (52% of the bursts) to approximately 60 s. 
For one burst the normalisation was found to be con- 
stant from 2-77 s after the start of the burst; the median 
duration was 48 s. Over the interval during which the 
blackbody normalisation was constant, the blackbody 
temperature typically decreased from 2.3 to 1.6 keV, 
while the flux decreased from w 3 x 10 -8 erg cm -2 s -1 



to 



5 x 10 9 erg cm 2 s 



5 The excluded bursts are #8, 15, 21, 
Table 5 of G08 



29, 39, 42, and 44 from 



The (Abb) values were not significantly correlated with 
either the duration over which the average was calcu- 
lated, nor the maximum time out to which the average 
was calculated. Interestingly, the earliest time to which 
the constant fits could be extended (« 1-3 s after the 
burst start) was within the «5s burst rise. For several 
of the bursts, the time at which the normalisation reaches 
the mean value corresponds approximately to a change 
in slope of the burst rise, seen in several of the bursts 
(e.g. Fig. [T]). This feature is similar to the "b ump" 
seen in model lightcurves (e.g. iHeger et al J 12007). Al- 
though in the simulations this feature must arise due to 
some variation in the rate of increase of thermonuclear 
energy production (as the 1-D models cannot account 
for lateral propagation), the coincidence of the observed 
change in slope with the normalisation achieving its mean 
value suggests that this point marks where the effects of 
spreading cease. Further analysis of the detailed shape 
of the burst rises may help to clarify this situation. 

We found significant variations in the mean normalisa- 
tion (Abb), both between and within different observa- 
tional epochs. For the bursts within each epoch (1997-8 
and 2000-7), the normalisation varied by 5 and 4%, re- 
spectively. This variation was highly significant; the % 2 
values for constant fits were 218 (for 6 degrees of freedom) 
and 1007 (for 59 degrees of freedom), respectively. In 
terms of the inferred radius, this implies systematic un- 
certainties of order 2-2.5%. Additionally, a much larger 
variation was measured between the two epochs. The 
first seven bursts, observed between 1997 November and 
1998 June, had (Abb) values substantially larger than 
the remaining bursts (observed from 2000 June onwards). 
The mean and standard deviation for the normalisations 
of the 1997-8 bursts was 118 ± 6 (km/diok pc ) 2 while for 
subsequent bursts was 101 ± 4 (km/ciiokpc) 2 - A two- 
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Figure 3. Mean blackbody normalisation values for 67 thermonu- 
clear bursts observed by RXTE from GS 1826—24. The mean and 
standard deviation for the two observation epochs (1997—8 and 
2000-7) are indicated. Note the marked discrepancy between the 
mean values for the two distributions; the K-S statistic indicate 
that they are discrepant at a significance of 1.3 X 10~ 6 , equivalent 
to 4.7a. 



sided Kolmogorov-Smirnov test confirms that these dis- 
tributions are discrepant at the 1.3 x 10 -6 significance 
lev el (equivalent to 4.7a) . This variation was also noted 
bv [Galloway ct al. (2004]), who reported instead a cor- 
relation between the persistent flux and the blackbody 
normalisation for a subset of the bursts analysed here. 

The observed variation is unlikely to arise from any 
instrumental effect, as the PCA is precisely calibrated 
to maintain stable flux measurements for calibration 
sources over the entire mission lifetime. As we discuss 
below, the discrepancy is also unlikely to result from vari- 
ations in the neutral column density n# as a function of 
epoch. Closer examination of the spectral variation in 
the two groups of bursts provides a possible explanation. 
In contrast with the example burst discussed in the pre- 
vious section, and other bursts observed in 2000-7, the 
constant fit interval for the bursts observed in 1997-8 be- 
gan later than in the 2000-7 bursts; typically 10 s after 
the burst start for the 1997-8 bursts, compared to 3 s 
after the burst start for the other bursts. To put this 
another way, there was additional variation in the black- 
body normalisation during the burst rise for the 1997-8 
bursts that prevented extension of the constant interval 
to the same point as in the later bursts. This is illus- 
trated in Fig. |4l which compares the blackbody normal- 
isation averaged over the 1997-8 bursts with that of the 
2000-7 bursts. Remarkably, however, the behaviour of 
the normalisation in the two groups of bursts is virtually 
identical during the burst rise; the discrepancy sets in 
from 10 s after the burst start, with the normalisation 
of the 1997-8 bursts gradually increasing over another 
ps 10 s to a higher level. 
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Figure 4. Comparison of averaged blackbody normalisation pro- 
files for bursts from GS 1826-24 measured in 1997-8 (#1-5 of 
G08) and 2000-7 (#9, 10, 11, 12, 13, 16, 17, 19, 20). The verti- 
cal dashed lines indicate the time of maximum flux for each set of 
bursts. Note the agreement in the normalisation throughout the 
burst rise and maximum, and the increasing discrepancy from 10 s 
after the burst start. The inset shows the corresponding variation 
of the averaged burst flux. 

3.2. KS 1731-26 

RXTE observations of KS 1731-26 in 2000 August- 
September detected 14 bursts, 8 in August and 6 in 
September (these are bursts #14-21 and 22-27 in G08, 
respectively). Although the RXTE/PCA observations 
were interrupted regularly due to the satellite orbit and 
other scheduled observations, the times of the detected 
bursts were consistent with a regular recurrence time. 
The shortest burst separation measured by RXTE dur- 
ing each month of data was « 2.5 hr; the longer sep- 
arations were consistent with integer multiples of this 
value, indicating regular bursts where intervening events 
were missed within data gaps. We measured the aver- 
age recurrence times based on linear fits to the burst 
arrival times measured by RXTE, as 2.577 ± 0.011 hr 
and 2.636 ± 0.003 hr for August and September, respec- 
tively. One further burst, observed on 1999 Aug 26 (#13 
in G08), exhibited a lightcurve consistent with the 14 
observed in 2000, and we included it in this sample. 

The regular bursts featured similar long (w 5 s) rises 
and decays (« 60 s) as those typically observed from 
GS 1826-24 (Fig. [5). The peak count rate was 
2300 count s -1 PCU -1 , so that deadtime corrections are 
comparable to those of GS 1826—24 (see However, 
the recurrence time for the bursts from KS 1731—26, 
at ss 2.6 hr, were significantly shorter than has been 
observed from GS 1826—24 in y ears of RXTE observa- 
tions (e.g. lThompson et al.1120081 ). Based on the assumed 
distances for the two sources (6 kpc for GS 1826—24 
and 7.2 kpc for KS 1731-26; see G08), the bursts from 
KS 1731—26 reached significantly higher luminosities of 
(1.8-1.9) x 10 38 ergs -1 , and the burst durations were also 
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Figure 5. Example burst observed by RXTE from KS 1731-26 
on 2000 September 29 14:08:35 UT (#24 in G08), illustrating the 
constancy of the blackbody normalisation over a significant extent 
of the burst rise and tail. Panel descriptions are as for Fig. [l] The 
reduced-x 2 = xt f° r the constant fit is 1.39, for v = 79 degrees of 
freedom. 

significantly shorter (r = 24 s compared with » 40 s for 
GS 1826-24). A shorter burst duration for KS 1731-26 
suggests a smaller fraction of hydrogen in the burst fuel, 
although the shorter recurrence time also allows less hy- 
drogen to be consumed by steady burning prior to the 
burst. 

As with GS 1826—24, we measured the best-fit normal- 
isation (-Kbb) from the time-resolved blackbody spectral 
fits over the longest possible time interval without ex- 
ceeding the 3er confidence limit. With the shorter bursts 
from KS 1731—26, the constant -Kbb fits extended typ- 
ically over the range 2-30 s after the burst start. Over 
this time interval, the blackbody temperature dropped 
(typically) from 2.5 to 1.8 keV, with the flux dropping 
from 3 x 10~ 8 crgcm _2 s _1 to 7 x 10~ 9 ergcm _2 s . 

The fitted (-Kbb) values of the long bursts from 
KS 1731—26 exhibited significant variability; a fit with 
a constant model gave a x 2 = 112.7 for 14 degrees of 
freedom, indicating variability at the 8.5cr level. The 
mean value was 80 ± 2 (km/d 10kpc ) 2 , with standard de- 
viation between the various measurements of 2.6%. The 
mean blackbody normalisation of the burst on 1999 Aug 
26, measured between 1.25-36 s after the burst start at 
79.7 ± 0.8 (km/<iiokpc) 2 , was fully consistent with the 
bursts observed in 2000 despite being observed a year 
earlier. 

In an independent analysis of the bursts observed from 
KS 1731-26 with RXT E (of which t he bu rsts we anal- 
yse here are a subset), IGuver et al.l ()201lD measured a 
mean blackbody normalisation of 88. 4±5.1 (km/diok pc ) 2 j 
and found weakly significant variation in the normali- 
sation as a function of the burst flux. Although this 
value is not significantly different from the normalisa- 
tion we derive from the subset of bursts analysed here, 
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Figure 6. Sensitivity of the blackbody normalisation in spectral 
fits of RXTE data to the assumed column density njj. The fits 
were calculated from a simulated spectrum with kT = 2.1 keV, 
including pre-burst persistent and instrumental background from 
burst #15 of KS 1731—26. The square symbols show the fitted 
value of the blackbody normalisation (left-hand y-axis) as a func- 
tion of the assumed column density n h ■ The lower half of the plot 
shows the corresponding \ 2 values; each fit is statistically accept- 
able. The choice of the nu value can clearly introduce a systematic 
bias to the normalisation values, of the order a few per cent. 

we suggest that the difference arises from a systematic 
effect due to the value adopted for the neutral column 
density. For the analysis presented here, we assumed 
riH = 1.3 x 10 22 cm~ 2 (from multi-epoch Chandra and 
XMM-Newton spec tra of the source during quiescence; 
iCackett et aLl l2006) while IGuver et al.l (|2011[) adopted a 
value of 2.98 x 10 22 cm~ 2 , the mean of best-fit values 
from fits to the burst spectra measured by RXTE. For a 
simulated blackbody spectrum with known kT and nor- 
malisation, adopting a realistic model for the pre-burst 
emission as background, the fitted value of the normal- 
isation depends linearly on the assumed nn (Fig. [6]). 
That is, over (under) estimating the nu value assumed 
for the fits will have the effect of over (under) estimat- 
ing the -Kbb- For the representative parameters chosen, 
the slope gives approximately 6.5 (km/(iiokpc) 2 for every 
additonal 10 22 cm~ 2 in column that is adopted. The dif- 
ference between the adopted values for the two analyses 
could account for an offset in the normalisations of up to 
11 (km/c?iokpc) 2 , sufficient to explain the discrepancy. 

4. DISCUSSION 

We found significant variations in the blackbody nor- 
malisation (-Kbb) averaged over tens of seconds of the 
burst tail in a homogeneous sample of regular, consistent 
bursts from GS 1826-24 and KS 1731-26, both within 
and between observation epochs. We found that the vari- 
ation within epochs was 4-5% (2.6%) for GS 1826-24 
(KS 1731—26), while the variation between epochs can be 
as large as 17% (for GS 1826—24). In terms of inferred ra- 
dius, this corresponds to variations of 2-2.5% (1.3%) and 
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8%. The va riation within each epoch is within the range 
reported bv lGiiver et al.l (|2011h in their more comprehen- 
sive study of bursts from several sources, although the 
8% variation between epochs observed for GS 1826—24 
is somewhat larger. We stress that this uncertainty is 
solely related to the measurement of the blackbody nor- 
malisation in bursts, and (as suggested by Zll) may be 
exceeded by other sources of uncertainty, for example the 
degree of anisotropy of burst emission. 

The 1997-8 bursts from GS 1826-24, which exhib- 
ited longer recurrence times and reached slightly higher 
fluxes, also exhibited larger mean normalisations (K^b). 
Comparison of the normalisation time-series averaged 
over subsamples of bursts from the two epochs show that 
the normalisations were essentially identical during the 
burst rise and through the burst peak, with deviations 
becoming apparent between 10-20 s after the burst start. 
The reproducibility of the -Kbb evolution throughout the 
burst rise and peak suggests that the system geometry, 
atmospheric composition and temperature (and hence 
f c ) were essentially identical over that interval, and what- 
ever physical condition gave rise to the elevated -fTbb in 
the 1997-8 bursts set in after the burst peak. We con- 
sider three possible mechanisms to give rise to the vari- 
ation. 

First, it may be that the amount of neutral mate- 
rial close to the neutron star was reduced, so that the 
neutral column density nu decreased during the 1997-8 
bursts, leading to an overestimation of iibb (see e.g. Fig. 
[6|). Simulations adopting the pre-burst emission from 
GS 1826-24 (as was done for KS 1731-26; see SgU) in- 
dicate that in order to overestimate the (K^) by the re- 
quired amount would necessitate decreasing the total uh 
value by 1.9 x 10 22 cm~ 2 . Since this value is almost five 
times the assumed line-of-sight value for GS 1826—24, 
we can rule out this mechanism. 

Second, it is possible that the degree of anisotropy 
of the burst emission changed in reponse to a varia- 
tion in the accretion geometry (e.g. lFuumotol[l988f) . per- 
haps triggered by the burst. Such a change might be 
expected to be reflected in the X-ray spectrum. The 
broadband (w 1-100 keV) spectral distribution of the 
persist ent emission in GS 182 6—24 has been well stud- 
ied by iThompson et ail (|2008l ). in the context of deter- 
mining the best possible estimate of the bolometric flux 
(and hence the accretion rate). Those authors found 
substantial changes in the spectrum with epoch, most 
notably during 2003 April, when evidence for an addi- 
tional soft (< 1 keV) component was found. However, 
the mean normalisation from the 6 bursts observed in 
2003 April was (102 ±4) (km/di k P c) 2 (see Fig. 0J), fully 
consistent with the mean from the other bursts from 
2000 onwards (but excluding the 2003 April bursts), of 
(101 ±4) (km/diokpc) 2 - Conversely, there is no evidence 
for a substantially differ ent X-ray spectrum du ring the 
1997-8 observations (see IThompson et al.ll200l Fig. 4) 
The lack of correspondence between the broadband spec- 
tral shape and the blackbody normalisation from the 
bursts seems to contraindicate an influence of the burst 
anisotropy on the normalisation, although cannot rule it 
out. On the other hand, the variations in the persistent 
X-ray spectrum might be transient, present only when 
the elevated iibb values were measured, beginning 10 s 



after the burst start. Such variations might be expected 
to lead to incorrect pre-burst emission subtraction later 
in the burst, although this is not observed consistently. 

Third, we consider the possibility that the f c value 
changes to produce the variation in the measured 
after the peak in the 1997-8 bursts from GS 1826-24. 
Modelling studies indicate that f c depends on fixed 
parameters such as the neutron star surface gravity, 
but also the composition and eff ective temperature 
of the scattering atmos phere (e.g. iMadej et al.l 120041 : 
ISuleimanov et all 1201 lbf ) . There is evidence in other 
sources that radius-expansion bursts can remove the 
outer, H-rich layers of the photosphere, leading to a 
change in the atmosph eric composition during the burst 
(Gall oway et al.l 120061 ); no such effects have been sug- 
gested from non-radius expansion bursts. Nevertheless, 
interpreting the maximum variation in (-Kbb) in the 
1997-8 bursts, compared to the mean for the 2000-7 
bursts, as a variation in f c implies a maximum vari- 
ation of 12% (8% in the mean). This result suggests 
that the assumption t hat f r is constant during b ursts is 
not always true. While ISuleimanov et al.l (|2011bl ) predict 
patterns of variation of f c with burst flux, the epoch-to- 
epoch differences in the burst lightcurves in GS 1826—24 
are relatively subtle, and do not appear sufficient to give 
rise to the inferred variation in f c . 

It is possible that the peak blackbody flux (or tempera- 
ture) serves as the discriminant which results in elevated 
blackbody normalisations later in the burst. The 1997- 
8 bursts reached maximum fluxes about 8% higher on 
average than for the 2000-7 bursts. The discrepancy in 
the maximum temperature reached was proportionately 
smaller. The samples of bursts studied here were deliber- 
ately selected for the consistency of their lightcurves and 
regularity of their recurrence times. For other samples 
of bursts, which are typically much more heterogeneous, 
systematic errors in measurements of the blackbody ra- 
dius of order > 8% should be assumed, unless the varia- 
tion in f c can be modelled. 

The mean blackbody normalisations measured during 
the later (2000-7) bursts from GS 1826-24, and the reg- 
ular bursts from KS 1731—26, were consistent with a con- 
stant value over several tens of seconds. This constancy 
was observed independent of the specific method for de- 
termining the mean value, although the choice of method 
can introduce small biases. In particular, comparisons of 
joint fits to the spectra with averages of the fitted nor- 
malisations show that in general the former method ar- 
rives at shorter durations for the constant normalisation, 
likely because the spectra are not (en masse) statistically 
consistent with blackbodies. Trends in the blackbody 
normalisation late in the burst tail, coupled with these 
marginally deviant spectra, likely give rise to systematic 
errors of a few percent between the two methods. The 
choice of time binning strategy does not have as large an 
effect. 

This constancy of the blackbody normalisation was 
maintained despite significant decreases in the black- 
body temperature, of (typically) 30%, and decreases in 
the flux of 60-75% over the same time interval. Such 
a lack of variation in the blackbody normalisation im- 
plies constraints on the relative degree of spectral dis- 
tortion over this temperature range, suggesting one of 
two situations. Either, any variation in in the spectral 
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distortion factor f c as a result of the varying effective 
temperature (as indicated by the decreasing blackbody 
temperature) is exactly balanced by some other varia- 
tion, e.g. a change in the emitting radius as the burst 
flux decreases; or, that the color temperature correction 
f c is also constant throughout the interval in which 
is constant. The former explanation is rather contrived, 
particularly considering that in the burst tail the burning 
front is expected to have already spread to the entire sur- 
face area of the neutron star, and no further increase (or 
decrease) in burning area is expected. Thus, these mea- 
surements suggest that for the range of effective temper- 
atures spanned in the burst tails, the color-temperature 
correction f c is roughly constant. This conclusion is dif- 
ficult to rec oncile with the prediction s of atmosphere 
models (e.g. ISuleimanov et all I2011bl ). which indicate 
significant variations in f c over most of the flux range 
spanned during the burst. The distance for G S 1826—24 
is thought to be « 6 kpc (Hc-g er et al.l 12007). at which 
an Eddington-limitcd burst would be expected to reach 
F/F Edd = 3.7 x 10~ 8 (6.3 x lO" 8 ) ergcm^s- 1 for an 
H-rich (pure He) atmosphere (e.g. G08). This would 
imply that the range of flux over which the burst nor- 
malisation is found to be constant is 0.16-0.75 (0.10- 
0.46) F/i^Edd- Th e greatest variation in th e normali- 
sation predicted by ISuleimanov et al.l (|2011aD is outside 
these ranges, although significant variations would yet 
be expected (particularly for the solar metallicity mod- 
els; see also Zll). 

The w 3-5% fractional variation in (Kbh) observed 
from GS 1826-24 and KS 1731-26 within each obser- 
vational epoch was smaller than that seen in two radius- 
expansion b ursts observed from EXO 1745—248, of 25% 
(|Qzel et al.l I2009T) . This observation perhaps suggests 
an explanation of the variation. EXO 1745—248 exhib- 
ited during it's 2000 outburst an initial period of strong 
variability, reminiscent of dip ping behaviour observ ed in 
high-inclination systems fe.g. iGallowav et al. 2008]) ■ No 
such variability has ever been observed from GS 1826—24 
or KS 1731—26, suggesting that perhaps the inclination 
is lower in those systems than in EXO 1745—248. If 
system inclination is the main factor in determining the 
variation in apparent blackbody normalisation, a possi- 
ble explanation is the reprocessing of some fraction of 
the burst flux off an accretion disk whose projected area 
varies with time. In 4U 1728—34, the ti mescale inferred 
for th e variation was several tens of days ([Galloway et al.l 
12003ft . much longer than the expected orbital period 
of this system (e.g. IGallowav et al.l feOlO'). However, in 
GS 1826-24 and KS 1731-26, the variation is on a much 
shorter timescale, comparable to the recurrence time of 
the bursts themselves (hours). For comparison, the dis- 
crepant bursts from EXO 1745—248 were separated by 
8.5 d. We note that such an explanation fails to ac- 
count for the factor of ~ 2 diff erence in normalisation i n 
the bursts from 4U 1724-307 (|Suleimanov et al.M2011ah . 
which does not show dips and thus is unlikely to be at 
high inclination. However, those bursts also exhibited 
markedly different timescales, indicative of a varying ac- 
cumulated fuel reservoir at ignition; we suggest instead 
that different physical conditions (temperature, compo- 
sition) gave rise to the difference in the measured black- 
body normalisation. 



The shorter timescale of variation in (ifbb) f° r 
GS 1826-24 and KS 1731-26 suggests that there may 
be an orbital component of the variation. The specific 
value of the blackbody normalisation depends upon the 
assumed neutral column density, as illustrated in Fig. [51 
Thus, an orbital variation in the line-of-sight column den- 
sity, perhaps arising from cool clouds of material above 
the point of contact of the accretion stream with the disk, 
will manifest as a variation in the blackbody normalisa- 
tion on the same timescale. In order to test this hy- 
pothesis, we calculated a Lomb-normalised periodogram 
on the blackbody normalisation measurements as a func- 
tion of time. Recall the orbital period in both sources is 
unknown, but for GS 1826—24 is thought to be around 
2 hr; based on the typical orbital periods for other burst 
sources, we searched a frequency range of 0.5-24 hr. For 
GS 1826—24, we divided the measurements from the two 
epochs through by the appropriate mean, and found a 
peak Lomb power of 10.48; for KS 1731—26, we found 
a peak Lomb power of 5.29. Neither of these detections 
is significant. For KS 1731—26, the known sensitivity of 
the normalisation measurement to discrepancies between 
the assumed and true value of nn indicates that a vari- 
ation of 0.3 x 10 22 cm~ 2 in Uh could account for the 
variation in the measured blackbody normalisation. For 
GS 1826—24, the measurements are slightly more sensi- 
tive to discrepancies in nu , so that a variation in njj of 
0.45 x 10 22 cm -2 could account for the variation in the 
blackbody normalisation within each epoch. If orbital or 
longer-timescale variations in nn are driving the varia- 
tion in the blackbody normalisation, it may be possible 
to verify through time-resolved high-spectral resolution 
measurements of absorption edges in the X-ray spectra. 
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